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ABSTRACT 


Stiffened plates are designed for minimum volume subject 
to constraints on: Von Mises maximum stresses, nodal dis- 
placements, height to thickness ratios of frame webs, and 
width to thickness ratios of frame flanges. Design variables 
are plate thicknesses and stiffener dimensions. 

A finite element analysis program is developed for the 
design of stiffened plates usiny numerical optimization 
technigues. The program may be used as a stand alone analy- 
Sis tool or may be coupled to an optimizer of user's choice. 
Rectangular plate elements and frame elements are used for 
the idealization of stiffened piates. 

Design examples are presented tc demonstrate the design 


method. 
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Lt. INTRODUCTION 


The development of high-speed digital computers has made 
possible significant changes in the structure design process. 
One of these is the availability of various mathematical 
programming methods for use in design optimization. The 
computer's speed allows the designer to now consider a much 
wider range of design alternatives. The optimization proce- 
dure provides a means of systematically choosing from among 
these alternatives based on some predetermined rational 
criterion. 

Even when the selected numerical method is able to arrive 
at the optimum design, tne result is only as good as the 
design model. Here there is an even greater need for design 
experience and sound engineering judgment. The design model 
must be carefully develcped to realistically represent the 
design in question. 

The finite element model for stiffened plate is shown in 
Fig. 1.1, where the plate has been descritized by rectangular 
plate elements and the stiffener by frame type elements. The 
eccentricity of stiffener is transformed to the linked nodal 
point by applying a linear equation that realtes displacement 
degrees of freedom. 

The purpose of this thesis is to develop a finite element 
analysis program for stiffened plates, and to design the 


optimum stiffened plate by coupling two programs; the 
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analysis program and ageneral purpose non-linear optimizer, 
COPES/CONMIN [Refs. 1-2]. 

The remaining chapters of this thesis are outlined as 
follows: 

Chapter II presents briefly the finite element method 
used in the analysis program for stiffened plates. 

Chapter III presents the basic concepts of the cptiri- 
zation methods used in the COPES/CONMIN. 

Chapter IV presents design examples. 


Chapter V offers conclusions and recommendations. 
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ie Fin teVEcEMeENT METHOD 


A. GENERAL 

The finite element method is now well established as an 
engineering tool of wide application. The fundamental 
concept of the finite element method is that any continuous 
quantity, such as pressure or «displacements, can be approxi- 
mated by a discrete modei composed of a set of piecewise 
continuous functions defined over a finite number of sub- 
domains. The piecewise continuous functions are defined 
using the values of the continuous quantity at a finite 
number of points in its domain. 

The formulation of the fin..te element method can be traced 
to energy procedures, principally the minimum potential energy 
method and the minimum complimentary energy method. The 
minimum potential energy method 1s associated with assumed 
displacement parameters as unknowns and is usually termed 
the "displacement" or "stiffness" method. On the other hand, 
the minimum complimentary energy method dealt with stress 
parameters and is termed the "flexibility" or "force" approach. 
The ease with which a continuous displacement pattern can be 
prescribed (compared to the alternative approach of forming 
an equilibrating internal force field) has aided the wide- 
spread use and development of the finite element displace- 
ment approach. The displacement model and the stiffness 


analysis are employed in the analysis program developed here. 


3 





This chapter will briefly present some of the general 
concepts of the finite element method used in the analysis 


program. 


B. FINITE ELEMENT DISPLACEMENT APPROACH 

The displacement formulation involves derivation of the 
stiffness matrix of each individual element. The stiffness 
matrix of the entire assembled structure is then obtained by 
Super-position. This matrix, along with the prescribed 
displacement boundary conditions and loads, is used for the 
solution of displacements and stresses. 

1. Element Analysis 

For the structural applications at least, the 

governing equilibrium equations can be obtained by mini- 
mizing the total potential energy of the system. The total 


potential energy, 1 can be expressed as 
il 
eo Vemma tedy =) ou “qds (ZBI) 


where d and € are the stress and strain vectors respectively, 
2 the displacements at any point, £ the body forces per unit 
volume and gq the applied surface tractions. Integrations are 
taken over the volume V of the structure and loaded surface 
eaea 5. 


Siew tiuse termeon tne paght hand side of Eq. (2-1) 


represents the internal strain energy and the second and 
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third terms are respectively the work contributions of the 


body forces and distributed surface loads. 


In the finite element displacement method, the basic 


steps for derivation of the element stiffness matrix are: 


a. 


Express the displacements to have unknown values only 
at the nodal points, so that the variation within any 
element is described in terms of the nodal values by 


means of interpolation functions. Thus 


.2.4) 


or 

lI 
ie 
ae 


where N is the set of interpolation functions termed 
the shape functions and ue is the vector of nodal 
displacements of the element. 

Express the strains within the element from the 


element nodal displacements as 


m 
lI 
2? td 
tC 
— 
nh 
ee) 


where B is the strain-displacement matrix generall 
composed of derivatives of the shape functions. 
Express the stresses relating to the strains by use of 


an elasticity matrix D containing the appropriate 


material properties, as follows: 


(2.4) 


Q 
II 

SEG, 
2m 


Establish the equilibrium equation of element. 


i) 





Provided that the element shape functions have been chosen 
so that no singularities exist in the integrands of the 
functional, the total potential energy of the continuum will 
be the sum of the energy contributions of the individual 


elements. Thus 


where T represents the total potential energy of element e 


mach, OM use of Eq. (2.1), can be written 
tT = a eee DD hia dy amg( er uN, LE dV 
V V 
2 2 
-~f ul yl g as (ee) 
Se ae ea 
e 


where is is the element volume and S, the loaded element 
surface area. Performance of the minimization for the 
element e with respect to the nodal displacements u\for the 


element results in 


art T e@ T AR 
meee {| 8(8 DE) ave= jf N £advV.-f N qds, 
Ju ee a Sa 
© Se S 
=k tice ~ aD 
where 
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is termed the element stiffness matrix. 
2. Direct Stiffness Method 

The real elastic structure is now represented by a 
finite number of small, discrete elements. Once their 
app:roximate behaviors, identified by their individual stiff- 
ness matrices k“of Eq. (2.9), have been established, the 
stiffness matrix K for the complete structure is obtained by 
the proper caeien of each element stiffness matrix in the 
structure. The summation of the terms in Eq. (2.7) over all 
the elements, when equated to zero, results ina system of 
equilibrium equations for the complete continuum. This 
assembly process is known as the Direct Stiffness Method. 
These equations are then solved by any standard technique to 
yield the nodal displacements. Note that K is symmetric and 


positive-definite. 


fe FINITE ELEMENTS USED 
The finite elements used in the analysis program will be 


described briefly in this section. 
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Mes Reeeenogubar Plate Element 


The rectangular plate element used here is illus- 
trated in Fig. 2.1, where each nodal point has 6 degrees of 
freedom, 3 transitional displacements and 3 rotational 
displacements. This element with 4 corner nodal points has 
the element stiffness matrix of order 24x24. The corres- 
ponding displacements of each node Uy rUsr--- Ung will be 
taken to be positive in positive directions of the xyz- 
coordinates. 

The main assumptions in the method are that displace- 
ments are small compared with plate thickness, the stress 
normal to the midsurface of the plate is negligible, and 
normals to the midsurface before deformation remain straight 
but not necessarily normal to the midsurface after deformation. 

The assumed displacement functions will be taken to 
be linearly varying in the plane of the element. These 
displacement functions will ensure both deflection and slope 
compatibility of the adjacent elements. The stiffness matrix 
of the rectangular plate element, which is matrix product 
B*DB of Eq. (2.9) integrated over its volume, is summarized 
by Przemieniecki [Ref. 3]. 

Note that element stiffness matrices are formed 
directly in the global coordinate system so that no trans- 
formations from local to global coordinates are required. 

2. Frame Eiement 
The frame element as a stiffener has 6 Gegrees of 


freedom for each node, such as those of plate elements. The 


dees: 
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basic assumptions are that the stiffener element is straight 
with uniform cross section capable of resisting axial forces, 
bending moments about the two principal axes in the plane of 
its cross-section, and twisting moments about centroidal axis, 
and that its deflection due to shearing strains are neglected. 
The stiffener element has 2 kinds of options which are illus- 
meaced in Fig. 2.3. The width and thickness of flange of 

the rectangular stiffener element may be referred as Zeros. 

In order to determine the stiffness property of a 
comp.ete structure, a global coordinate must be established 
for all unassembled structural elements so that all the 
displacements and their corresponding forces will be referred 
to this system. Since the element stiffness matrices K are 
imitlally calculated in local coordinates, suitably oriented 
to minimize the computing effort, it is necessary to intro- 
duce transformation matrices changing the frame of reference 
from a local to a global coordinate system. The first step 
in deriving such a transformation is to obtain a matrix rela- 
tionship between the element displacements u in the local 
system and the element displacements u’ in the global system. 


Te 


This relationship is expressed by the matrix equation 
16 Ie eee ta 6 Beg Zen) 


where T is a matrix of coefficients obtained by resolving 


~ 
~~ 


global displacements in the directions of local coordinates. 


The transformation matrix T 1S given by 


Ld 
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Pirgune 2.2 Frame Element 


zal: 


a. Rectangular type frame stiffener 





b. T type frame stiffener 





Figure 2.3 Stiffener Options 
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A 0 0 0 
0 A 0 0 
T = bs 
. 0 0 A 0 
0 0 0 A 
where 
g m n 
Ox OX OX 
A = Q m Ty 
~ Oy OY Oy 
AGE Moz Noz 


represents matrices of direction cosines for local x, y and 
z directions, respectively, measured in global system x, y 
and z. Similarly the relationship for forces is expressed 


oy the: matrix equation 
p = Tp (eae) 


wnere p is the force vector in the local system and p' is 


~w ~- 


the force vector in the global system. Matrix A is orthogonal; 


tnat is, r = r. Therefore 

Hoe : u ae) 
and 

p! = T's (Q22) 


23 


Let u' and u be two ways to describe the same virtual dis- 


~ 


placement. Virtual work is 


jl Ge Saree sul?" p (2.14) 
SO 
su(p'-T’p) = 0, pt = Tp (2.15) 
To transform the stiffness matrix, we start with on = p' 


and substitute from the preceding equations. 


oD tay ee ee ee ae 


Bamce this relation is presumed valid for any u‘, we conclude 


that the required stiffness transformation is 
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Pee OPFSETTING OF RIGID LINKS 

One of the most important advantages of the finite ele- 
ment technique is that an assembly of different structure 
elements such as plates and frames can be dealt within a 
Single coordinate analysis. Usually the neutral surfaces of 
plate and stiffener are not coincident: the stiffener 1s on 
one side of the plate. A standard preliminary treatment is 
to connect adjacent plate and stiffener nodes by a rigid 


link, so that degrees of freedom of the stiffener are replaced 


24 





by degrees of freedom of the plate. The usual assembly is 
then possible. The necessary transformation is now described. 
The stiffener element in Fig. 2.4 has the usual 12 degrees 
of freedom--6 at node A and 6 at node —E. With reference to 
these degrees of freedom, element load and stiffness matrices 
are Pp and k. Similar degrees of freedom are used at nodes l 
and 2 of plate element of rigid links A-l and B-2. The 
"master" degrees of freedom at node l end "Slave" degrees 


of freedom at node A have the relation 


Ur = Yi as) 
4A2 maze 
7A3 U3 | 
“a4 “14 : 
"AS "15 | 
YA6 "16 | 
where 
A, = i 0 0 0 te -2oy 
0 iM 0 ey G al 
0 0 i Pel ey 0 
0 0 0 Ll 0 0 
0 0 0 0 i 0 
0 0 0 0 0 1 
and hy = %) 7 Xr “ohh = ie mee and oes = Z1- Zy> 


ZD 


plate element 


frame element 


Figure 2.4 Rigid Links 
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A similar expression 1S written for link B-2 by replacing 
subscripts A and 1 by B and 2. The transformation vectors 
p and matrix k, associated with degrees of freedom at nodes 1 


st 


and 2, are 


p = Tp (2.19) 
k = rT Kt (P20) 
where 
entero es eed 
2 


The foregoing transformation makes the translational 
displacements depend on the rotational displacements and so 
introduces an unwanted quadratic field into the translational 
displacements. It is found that a frame-stiffened plate is 
overly flexible in a coarse mesh. Mesh refinement helps: 
error is reduced by a factor of four if the number of elements 
is doubled. The error can be eliminated through the addition 
of one more degree of freedom. Further details of the error 


reduction are described in [{Ref. 4]. 


2,9) 


Le Ober ZAT LON Lenn reues 


A. INTRODUCTION 

In this chapter some of the fundamental ideas and formu- 
lation methods of mathematical programming are introduced 
to understand the optimization method and the optimizer 
COPES/CONMIN. DeSign parameters used here are then dis- 
cussed. Fox [Ref. 5] and Himmelblau [Ref. 6] provide an 
extensive discussion of numerical optimization techniques 


and their application to engineering design. 


See DEFINITIONS 

In discussing the optimization methods, the following 
definitions will be useful: 

Design variables--the design variables are the numerical 
parameters for which values are to be chosen in producing a 
design. In a structural problem, they might be plate thick- 
ness, frame dimensions, etc. 

Objective function--the objective function, is the single 
valued function with respect to which the design is optimized. 
In a structural design problem, it might be the weight, volume 
or fabricated cost of the structure. The selection of an 
objective function can be one of the most important decisions 
in the optimum design process. 

Constraints--practical design problems are usually sub- 


ject to a series of constraints which must be satisfied in 
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order to produce an acceptable design. These constraints 
may be linear or non-linear. If a parameter is beyond the 
value of a specified value, it is said to be violated. 
Side ccnstraint--the side constraint is a constraint 
which restricts the specified range of a design variable 
for reasons other than the direct consideration of performance. 
Feasible design--the feasible design (or acceptable 
design) 1s a design which satisfies the specified constraints. 
Infeasible design--the infeasible design (or unaacep- 
table design) is a design which does not satisfy the speci- 
fied constraints. 
The general nonlinear optimization problem is expressed 


mathematically as follows: 


Minimize F(x) Gael) 


Subject to: 


Ses, <8 £Ob =. ..5m (32) 
x” < x. < xo fOr Ae = iy (32250 
i- i- 3. 


where the vector, x, is the vector of n design variables. 

The objective function, F(x), given by Eq. (3.1), as well 
lemene Constraint functions given by Eq. (3.2), may be linear 
or non-linear functions of the design variables. They may be 
eeooiacit OF implicit functions of x, but must have continuous 
first derivatives. If it is desired to maximize F(x), the 


“es 


minimization of the objective function is used since maximum 
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Di) can be treated minimum of -F(x). The function g.. (x) 
is the set of inequality constraints to be met. The m in- 


equality constraints must be satisfied to be a feasible 


design. Side constraints, ss and Kay are the lower and upper 
limitations placed on the design variables. Side constraints 


could be included in Eq. (3.2), but are treated separately 
for efficiency. 

The n-dimensional space spanned by the design variables 
Xs is referred to as the design space. As stated previously, 
any design which satisfies the inequalities of Eq. (3.2) is 
referred to as a feasible design. If the design violates one 
more of the inequalities, it is said to be infeasible. The 
minimum feasible design is said to be optimal. 

Most nonlinear optimization programs update the vector 


of design variables by the iterative relationship 
= eee esd (3.4) 


where gq is the iterative number, vector S is the direction 
of search in the design space, and a* is a scalar which 
defines the distance of travel in the direction S during the 
qth iteration. An initial design defined by xO must be 
defined and may be a feasible or infeasible design. 

The optimization process then proceeds in two steps. The 
first is the finding of S which will improve the design 
without violating constraints and the second is the determina- 


tion of a* which will improve the design as much as possible 
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when moving in this direction. The process is then repeated 
until there is no further design improvement:, indicating 


that this is the optimum attainable design. 


C. COPES/CONMIN 

The COPES/CONMIN optimization program is a general pur- 
pose, non-linear optimizer capable of handling large, con- 
strained problems. This includes the conjugate direction 
method of Fletcher and Reeves [Ref. 7] for ‘sanconstrained 
function minimization and a modification of Zoutendijk's 
Method [Ref. 8] of Feasible Directions for constrained func- 
tion minimization. It has been successfully used in connec- 
mem with structural optimization [Ref. 9], airfoil design 
(Ref. 10], aircraft synthesis [Ref. i1!, and numerous other 
engineering applications. 

It was necessary to develop a subroutine, ANALIZ, which 
for a given design, would analyze a stiffened plate, and which 
would be suitable for coupling with the optimizer. The 
common block GLOBCM is required to couple the analysis sub- 
memeine directly to the COPES/CONMIN. All variables, which 
are used as objective function, constraints and design varia- 
bles, must be listed in the common statement and the statement 
must appear in each subroutine the variables are used in. 

It is used by the optimizer as a catalog to identify where 
the design variables, objective function and constraints are, 


and what purpose they fulfill. 
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In order to execute the COPES/CONMIN program it is 
necessary to provide formatted data for COPES/CONMIN, followed 
by data fcr the ANALIZ program. Further details for the 


COPES data are explained in [Ref. l]. 


D. DESIGN PAFAMETERS 
Tne design variables, constraints and objective function 
used in the design process are discussed in this section. 
i) Pesigmevaintables 
The design variables are plate thicknesses, ti, and 
stiffener dimensions, being the web heights ha, the web 


thicknesses t the flange widths Wie and the flange thick- 


Wi’ 


nesses tees 1 = l1,ns, where ns is the number of the different 
stiffener dimension sets. Thus, the total number of design 
variables depends on the sets of plate thicknesses and the 
sets of cifferent stiffener dimensions. 

MeO SIeCe Ve Ch uncerOn 


Total structure volume is considered as the objec- 


tive function in the design process. 


™m ns 
OR ote 
i 


Volume = 
where a is the volume of the ith plate element, m is the 
number of plate elements, V 54 is the volume of the jth stiffener 
set, and N, is the number of stiffeners in this. The number 


@®= Staffener sets is ns. 


2 


——<¢ 





Se —COlistraints 
Design constraints are Von Mises maximum stress, 
nodal displacements, height to thickness ratios of frame 


webs, and width to thickness ratios of frame flanges. 


Stress: 
g@ JO - 1 < 9 ior ne 
iy Meas = ide sae 
S eae 7) 
poe - i < 90 = iene 
a) ier x — / 
J 1,nj 


where ae and ee are respectively the upper and lower sur- 
face stress at the node j of the element i, and one is the 
maximum allowable stress. The number of elements 1s ne and 
nj 1s the number of nodal points. 

Displacements: 


Os a ae =e Le, i 7 Eas 


where re is the displacement at node i in the coordinate 


Sreection j and u is the maximum allowable displacement 


max 
in coordinate direction j. The number of coordinates 1s 
ne. 


Height to thickness ratio of frame webs: 


i = 1,ns 


U1 
A 
>) 


hi/twi 7 = 


A 
© 


Te 


33 





where A, and) t are respectively the wek height and web 


Wi 
thickness of ith stiffener set. 


Width to thickness ratios of flanges: 


Wi/tes =o 95 po i= i,ns 


A 
en) 


Ih SE Oc 


where Ws and tes are respectively the width and the thickness 


of flange of the ith stiffener set. 
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IV. DESIGN EXAMPLES 


This chapter will present design examples of unstiffened 
and stiffened plates. All calculations were carried out on 
the IBM 3033 model 370 computer at the Naval Postgraduate 
School. Material constants used for all examples are Young's 
modulus, E = 2) aI psi, and Poisson's ratio, v = 0.3. Due 
to symmetry only a quarter of the plate was modelled in each 


example. 


A. DESIGN CASE I-~-THE UNSTIFFENED PLATE 

A quarter of the unstiffened plate simply supported, as 
shown in Fig. 4.1, was modelled using a 4x4 mesh of the 
plate elements. There are 10 design variables as illustrated 
in Fig. 4.2. The volume of the plate was minimized subject 
to a total of 278 constraints being stress and nodal defiec~ 
Etom Constraints, with Daas = 20,000 psi and Gens TP 2s 

im Case TA with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at the 

center of the whole plate. The final optimization results 
show that there are 3 critical deflection constraints o£ 
the nodes 3, 4 and 10, and that the plate thicknesses of the 
diagonal elements are much thicker than those of off-diagonal 
elements. The final volume is 363.916 in?, The results of 


this case are summarized in Table I for comparison with the 


other cases. 


5), 
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& @ ° 
(eae oe 
y t : plate thickness 
@ <= element 


@ : node 


Eugqure 45.) Design Case I--The Unstiffened Plate Model 
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bagure <4. 2 Design Case I--Design Variables 
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TABLE I 


Summary of Design Case I 


Variable Concentrated Distributed 
Load Load 
m1) the SO 02552 
E{2) IL 10) il O263 
ies) OPS 78 0.589 
(4 ) bso) 02729 
ie ) OAL AL 0.086 
6 ) OZ 0 OS SP cis, 
te 7 ) ORLA 0.086 
ito) O21 20 U0 
tee 9 ) On 120 0.086 
iLO) O20 ORCS 5 
Vol Bios o 6 70 2 os0 


t: plate thickness (in) 
Vol: volume (in”) 


fee case TB ith Unt torme Distributed Load 

Pe wil rOorml Gl strrbuted GlOad ja 0.276 lbs/in* was applied 
over the whole plate. The final optimization results show 
that there are 6 critical deflection constraints of the 
nodes 3, 4, 5, 9, 10 and 15, and that the plate thickness 
has a similar trend to case IA where diagonal elements are 
thicker than those of the off-diagonal elements. The final 
3 


volume is 170.396 in The results are summarized in Table I 


with those of case IA. 


Error olGN CASE TI=--THE STIFFENED PLATE 

A quarter of the stiffened plate simply supported, as 
shown in Fig. 4.3, was modelled 4x4 mesh of the plate ele- 
ments and 12 stiffener elements. 


in i @ |} 
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Figure 4.3 Design Case II--The Stiffened Plate Model 
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iPemcacce lia Receangular Iwoe Frame Stittener 


There are 17 design variables: a plate thickness, 
and 8 sets of stiffener heights and thicknesses, as tabu- 
lated in Table II. The volume of the stiffened plate was 
minimized subject to a total of 222 constraints: stress, 
stiffener height to thickness ratio, and nodai deflection 
constraints. 

a. Case IIA1l with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at 
the center of the whole plate. There are 14 crit:.cal con- 
Straints: the height to thickness ratios, and 6 deflections 
of the nodes 4, 5, 10, 29, 30 and 31. The final volume is 


M0017 in> 


Wht fa Undrorm plate-thackness of 0.051 i1n,.° The 
results are summarized in Table III for comparing with the 
other cases. 

b. Case IIA2 with Uniform Distributed Load 

A uniform distributed load, 0.278 Lb/in? was 

applied over the whole plate. There are 12 critical con- 
straints: the height to thickness ratios, and 4 deflections 
of the plate nodes 7, 8, 14 and 19. The final volume is 
65.984 in? with a uniform plate thickness of 9.059 in. The 
results of this case are summarized in Table III with those 
of the case IIA1L. The volumes of the above two cases are 


significantly reduced in comparison to those of the unstif- 


fened plate. 
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TABLE II 


Design Case IJA--Design Variables and Nodal Connectivity 


element design node number 

number variables Nl N2 N3 N4 
pl c J 6 2 
p2 te ja a 3 
joe. 2 3 8 4 
p4 c 4 9 10 S) 
p5 t 6 ba 12 i 
p6é ie Fi d2 ie 8 
p7 tz 8 ales 14 9 
ps8 is 9 14 LS 10 
ON, te iat 16 Bs Ik2 
plo t 2 ee 18 dae} 
> 1 ( 1s 18 We 14 
0) t 14 Ie, 2) ts) 
joes: e 16 pat a2 iy 
p14 t ay, ae 23 18 
PLS t im: 2S 24 Je, 
pl6é iS 18, 24 vis) 20 
sl t (1) ,n(l) 30 35 
s2 t,.(2) ,h (2) 35 36 
$3 t..(3) ,n(3) 36 CF 
s4 ty (4) ,h(4) oy, 38 
s5 t,. (5) ,h(5) 29 a) 
s6 t,,(6) ,n(6) 23 29 
s7 ty (7) ,n(7) 27 28 
s8 ty (8) ,h(8) 26 27 
s9 ty (5) ,A(5) 30 BL 
sl10 t,, (6) ,n(6) Suh S2 
sll ty (7) a7) Be 38 
sl2 t (8) ,n(8) 38 34 
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tees LET 


Summary of Design Case IIA 


variable concentrated adzstributed 
load load 
S Cpe sal CaO sy 2, 
tt) 0.438 O22 03 
ty (2) Os Sills: 0-168 
ty (3) OF 553 C2267 
th (4) Oro > 0.344 
ty (9) OF 0357 OL LS8 
t,, (6) O06 0.164 
ty (7) O02 5 Oly 1! 
ty (8) OOS Opa Sys) 
men ( 1.) 25189 1.044 
fet 2) 12/568 0.832 
ne 3 ) Pret sey 23:29 
h (4) 3.343 ee) O19 
rin ) 0.486 Omos? 
i 0) 0.482 Ong 2a 
ee) ) 0.478 O85 
macs) Oe477 O27 76 
Vol EO2 017 65.984 


plate thickness (in) 


frame thickness (in) 


pa ct 


frame height (in) 


Vol: volume Gn) 
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eee aseel in.) TVpcubeuame str tfener 


There are 33 design variables: a plate thickness, 8 
sets of the web heights and thicknesses, and 8 sets of the 
flange widths and thicknesses, as tabulated in Table IV. 

The volume of the stiffened plate was minimized subject to 
a total of 230 constraints: stress, web height to thickness 
ratio, flange width to thickness ratio, and nodal deflection 
Sonstraints. 
a. Case IIB1 with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at 
the center of the whole plate. There are 1Z critical con- 
straints: the web height to thickness ratios, the 3rd and 
8th flange width to thickness ratios, and deflection con- 
Straints at the node 5 of the element p5 and node 30 of the 
element sl. The final volume is 55.091 in? With a Undiorm 
plate thickness of 0.023 in. The results of this case are 
summarized in Table V for comparison with the other cases. 

b. Case IIB2 with Uniform Distributed Load 

PeUnizewrm distributed load, 0.278 ibs /in* was 
applied over the whole plate. There are 16 critical con- 
straints: the web height to thickness ratios except the 
8th one, the 3rd through 8th flange width to thickness ratios 
except the 6th one, and 4 deflections of the plate nodes 7, 


3 With a -UniLrorat 


8, 14 and 19. The final volume is 64.414 in 
plate thickness of 0.060 in. This volume is greater than the 
volume of the above concentrated load case, being different 


from the other cases which the latter is greater than the 
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TABLE IV 


Design Case II5--Design Variables and Nodal Connectivity 


element design node number 

number variables Nl N2 N3 N4 
onk 2 af 6 2 
joy 12 Z 7 8 S) 
ions tf 3 8 4 
p4 1 4 3) 10 5 
p5 ie 6 glee ie2 y 
pé 2 yi lee cs 8 
py 12 8 3 14 g 
ps8 ic 9 14 5 10 
pg t suse 16 any 2 
“ollie, c ie Bae, 18 i3 
en ie as 18 a2 14 
p12 te 14 aS 20 as 
Bi 3 1 16 2 22 ue 
p14 ie 7, Ze 23 18 
D5 ‘= 18 a3 24 19 
plé e ARS, 24 25 pa) 
sl ty(l),A(l),t,. (1) ,w(l) S0 Se 
S2 ty(2),h(2),t,.(1) ,w(l) 35 36 
$3 t.,(3) 2 (3) ,t¢ (3) ,w(3) 36 oy 
s4 ty (4) ,n(4),t, (4) ,w(4) 37 38 
s5 ty (5) ,n(5),t¢ (5) ,w(5) PS, 30 
sé ty (6),h(6),t,(6) ,w(6) 28 PL, 
s7 ty (7) ,n(7), te (7) ,wl7) 27 28 
s8 ty (8) ,n(8) ,t.(8) ,w(8) 26 ZT 
s9 ty(5),h(5),t¢.(5),w(5) 30 31 
s10 ty (6),h(6),t,(6) ,w(6) SA 32 
alph ty (7) ,n (7), te (7) ,w(7) Se a3 
sl2 ty (8) ,n(8),t.(8) ,w(8) a3 34 
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PAB ay 


Summary of Design Case IIB 


variable concentrated distributed 
load load 
te OOS) 0.060 
ty (1) Oris 3 0.149 
ty (2) ORO) 0.126 
ty (3) Gres / Ce 7 
ty (4) OZ 0.242 
ty (5) Oe ereil CeshO2 
ty, (6) O355 Oka 
ty (7) 0.316 Ores 
ty (8) 0.274 Ort 2 
in (1) ORE SS, 05.736 
lage 2 ) 0.386 O76 31 
ia (3) 0.680 On 76 
h (4) C305 a ey. 080 
ee 2 ) e902 0.509 
h (6) Lao’ 02566 
ee 7 ) due) 075629 
h (8) Peses Ono 02 
t, (1) OPO) 28k OiZoZ 
(2) 0.028 0.183 
t, (3) 6) OBE 0.308 
t, (4) 02026 0.460 
te (5) OS 76 C209 
t, (6) 0.799 Oe 23 
te (7) Ose alk On 257 
8) 0.460 0.201 
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TABLE V (Cont'd) 


variable concentrated dgseributed 
load load 
w (lL) ORCS O24 s 
w (2) 0.054 Ome oy, 
w (3) Oa057 02307 
w (4) 0.074 0.460 
w (5) OmS9 1 Or 2260 
w (6) 0.814 OnZ 35 
w (7) 0.685 22213 
w (8) 0.464 OF Zonk 
vol ao. 091 o4.4].4 
t: place thickness (in) 
tig: web thickness (1n) 
h: web height (in) 
te? flange thickness (in) 
w: flange width (in) 


Vol: volume (in?) 
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former. This suggests that the stiffener type and arrange- 
ment used here is particularly efficient for the concentrated 
load case. The results of this case are summarized in Table 
V with those of the case IIBl. The volumes of the above 
cases are Siqnificantly reduced with comparison to those of 
unstiffened plates. The volume of T type stiffener case 

is slightly smaller than that of rectangular case in the 
distributed load. On the other hand, the volume of T type 
stiffener case is much smaller than the rectangular case 


for the concentrated load. 


meno IGN CASE LIIl=--THE HATCH COVER 
A guarter of a hatch cover was modelled 4x4 mesh of the 
plate elements and 24 stiffener elements, as shown in Fig. 
4.4. 
ee ace el AST Rec bangulace ype: Prame Stiftener 
There are @ design variables: a plate thickness, 
and 4 sets of stiffener heights and thicknesses, as tabulated 
in Table VI. The volume of the hatch cover was minimized 
subject to a total of 274 constraints: stress, height to 
thickness ratio, and nodal deflection constraints. 
a. Case IIIA1 with A Concentrated Loac 
A concentrated load of 1,000 lbs was applied at 
the center of the whole plate. There are 4 critical con- 
straints: the 2nd and 4th height to thickness ratios, and 
the upper and lower stresses of the node 5 of the element 
p4. The dimensions of the 2nd stiffener set are negligibly 


small with comparison to those of the other stiffener sets. 
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CASE A 


O : plate element 
QO : stiffener element 
e : plate/stiffener node 





CASE 8B 





Figure 4.4 


SECTION A-A AND SECTION B-8 


Design Case III--The Hatch Cover Model 
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TABLE VI 


Design Case TIIA--Design Variables and Nodal Connectivity 


element 
number 


pl 
p2 
oe 
p4 
ope. 
po 
1 / 
ps 
py 
plod 
pill 
pl.2 
pl3 
p14 
B15 
plo 
sl 
S2 
$3 
s4 
s5 
$6 
s/ 
$8 
s9 
s10 
sll 
Suz 


design 
variables 


ROR CRueCae (er tu nr (ne tae ne tar tae a tr 


2 
1/2t,.(3) ,h(3) 
1/2t,(3) ,h(3) 
ty (1) ,h(l) 
t(1),h(1) 
ty (4) ,h(4) 
ty (2) ,n(2) 
ty (1) ,n(1) 
ty (4) ,n(4) 
ty(2),h(2) 
ty (4) ,A(1) 
ty (3) ,h (3) 
ty (3) ,h(3) 
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rr Oo moO DBD FP WD N FF 


N2 
6 
7 
8 
, 
bal 
12 
LS 
14 
16 
ee 
18 


node number 


N3 


Oo CG ~I 


Z 
rs 


hem WO O yD MO SF W N 





element 
number 
s1l3 
sl14 
ee 5 
sl6 
sl7 
s18 
sl19 
$20 
s2l 
S22 
S23 
$24 


TABLE VI 


d2sign. 
variables 
t(2) ,h(2) 
t..(2),B(2) 
ty (4) ,B(4) 
t.7 (3) ,h(3) 
1/2t,,(3) ,h(3) 
ty (4) ,n(4) 
t.7(3) ,f (3) 
1/2t,,(3) ,h(3) 
ty (4) hh (4) 
ty (4) ,n(4) 
ty (4), (4) 
ty (4) ,h(4) 


(COM “d) 
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N1 


36 
5 / 
34 
36 
38 
Shs, 
40 
41 
42 
43 
44 
45 


node number 


N2 


Cy 
39 
oy, 
40 
41 
42 
44 
46 
43 
44 
45 
46 


N3 


N4 





This suggests that the lst stiffener set is stiff enough to 
compensate the elements attached to the 2nd stiffener set. 


Meemeinal volume is 3230639 in- 


with a uniform plate thick- 
ness of 0.359 in. The results of this case are summarized 


maka ble VIL. 


TABLE VII 


Summary of Design Case IIIA 


variable concentrated distributed 
load load 
fs 0.4359 O7059 
ty (Ll) Orro lL 1 0.863 
ty (2) Ord 01 02002 
ty (3) 0.084 Osos 
ty (4) Ov OS 2009 
Fee GL} Lease a 750 
ha (2) 0.006 02001 
ne £5) ORS 1.484 
h (4) Ci0s O20? 
Vou B22 650 Bi esi4 


t: plate thickness (in) 


tig: frame thickness (in) 
h: frame height (in) 
Vol: volume (in?) 


b. Case IIIA2 with Uniform Distributed Load 
Amun form caisertourea load, 0.278 ey ah was 
applied over the whole plate. There are 6 critical constraints: 


the 2nd height to thickness ratio, and 5 deflections of the 


ol 





plate nodes 8, 9 and 14 and of stiffener nodes 32 and 37. 
With the same suggestion as the concentrated load case the 
dimensions of the 2nd stiffener set are negligibly smeil. 
The final volume is 87.314 in? with a uniform plate thickness 
of 0.059 in. The design results of the above both cases 
are summarized in Table VII with those of case IiIAl. The 
volumes of the above two cases are significantly reduced 
with comparison to those of unstiffened plates. 
eee ea cCePtrE. =. byeer! Frame Stilrrener 

There are 17 design variables: a plate thickness, 4 
sets of the web heights and thicknesses, and 4 sets of the 
flange widths and thicknesses, as tabulated in Table VIII. 
The volume of the hatch cover was minimized subject to a total 
of 278 constraints: stress, web height to thickness ratic, 
flange width to thickness ratio, and nodal deflection 
constraints. 

a. Case IIIB1l with A Concentrated Load 

A concentrated load of 1,000 lbs was applied at 

the center of the whole plate. There are 3 critical con- 
straints: the 3rd flange width to thickness ratios, and 2 
deflections of the plate node 5 and stiffener node 30. 
Negligibly small are the whole dimensions of the 2nd stiffener 
set and the flange dimensions of the 3rd stiffener set with 
comparison to the other dimensions. The final volume is 


a, 30 3 in? WEEN ea Und tOrmMe plate sEhickness, 0.083 in. 
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TABLE VILE 


Design Case IIIB--Design Variables and Nodal Connectivity 


element design node number 
number variables N1 N2 N3 N4 
pl te il 6 2 
p2 ic 2 7 3) 
oe) te 3 8 4 
p4 te 4 2) 10 5 
pS = 6 ile EZ ° 
po te zi lige ae: 8 
jou te 8 LS 14 S, 
ps te S) 14 ali 10 
p9 1 Jae 16 17 eZ 
p10 — 1 i 18 13 
jellal te 1a 18 io 14 
pl2 te 14 19 20 nS 
oulgs) b 16 aall 22 iy 
p14 te lew ie 23 18 
pl5 te is 25 24 IE) 
pl6 c Veg) 24 25 20 
sl 1/2t.(3) ,h(3),t,(3) ,1/2w(3) 26 27 

$2 1/2t..(3) ,h(3) ,t¢(3),1/2w(3) 20, 28 

s3 ty(1),n(1),t,(1),w(l) 28 Pu, 

s4 ty(1),h(1),t,(1) ,w(l) Zo 30 

$5 ty (4) ,n(4) ,t,(4) ,w(4) 26 Sil 

S6 ty (2), Rn(2),t_(2) ,w(2) 28 a2 

s7 t. (1),h(1),t,_¢(1) ,w(l) 30 5) 8 

s8 ty(4),hn(4),t<.(4) ,w(4) ul 34 

s9 ty(2),h(2),t,(2),w(2) 32 36 

sl10 ty(1),h(1),t<(1),w(l) es 38 

aie ti. (3),h(3),t— (3) ,w(3) 34 a5 

SZ ty (3) ,h(3) ,t¢(3) ,w(3) 33 36 


oye 


element 

number 
s13 
s14 
S15 
sl16 
sl17 
sl18 
sl19 
$20 
s21 
S22 
$23 
s24 


Eiepoevril (Cont'd) 


design 
variables 

Era IN 2) ne 102), wit 2) 

try (4) (4) bt, (4) ,w(4) 

eee eae) ren), wi) 
1/2t,,(3) ,-h(3),t¢(3) ,1/2w(3) 

ty (4) H(4) te (4) ,w(4) 

Sag ) yl) gale) pee 
1/2t. (3) 7n(3) ,t_ (3) pe 2 \ |) 

ty (4) -n(4) ,t,¢(4) ,w(4) 

ty (4) (4) ,t¢(4) -w(4) 

ty (4),H(4),t-(4) wld) 

try (4) (4) te (4) ,w(4) 
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Nl 


36 
57 
34 
36 
38 
oy, 
40 


42 
43 


45 


node number 


N2 


Si 
38 
By) 
40 
41 
42 
44 
46 
43 
44 
45 
46 


N3 


N4 





Dowcease llIp2 with Uniform Distributed Load 

A uniform distributed load, 0.278 lbs/in* was 
applied over the whole plate. There are 9 critical con- 
straints: the 3rd web height to thickness ratio, the 3rd 
flange width to thickness ratio, 4 deflections of the plate 
modes 8, 9, 13 and 14, and 3 deflections of the stiffener 
nodes 32, 36 and 37. Negligibly small are the whole dimen- 
Sions of the 2nd stiffener set, and the flange dimensicns o: 
the 3rd and 4th stiffener sets with comparison to the other 
dimensions. This suggests that the lst stiffener set is 
stiff enough to compensate the elements attached to the 2nd 
stiffener set. Additionally it can be suggested that the 
3rd and 4th stiffener sets are stiff enough to support the 
geezer l1oad with only the rectangular type ones. The final 
volume is 76.436 in? with a uniform plate thickness, 0.062 
in. The results of the above both cases are summarized in 
Table IX. The volumes of the above two cases are signifi- 
cantly reduced with comparison to those of the unstiffened 
plates and the volume of T type stiffener case is much 


smaller than that of the rectangular type case. 


ape 


TABLE IX 


Summary of Design Case IIIB 


variable concentrated Grstributced 
load load 
t 0.083 02062 
t (1) 0.284 0.284 
ty (2) OF Gs Oe 00s 
t (3) 0325 Oknecule 
ty (4) O02 p 02020 
lau GL, ee Oli 1.200 
Iee(-Z ) O01 8 0.020 
rie ) Lega S 3 1.430 
h (4) 0.1160 0.060 
t.(1) in0Od 1.000 
t- (2) 0-008 O2003 
t (3) 0.100 O00 
t.(4) 0.010 0.010 
w (1) 5 Se Ce 5 ere 8 
w (2) O1..01-0 0.010 
w (3) 02002 2003 
w (4) 0.930 0.001 
Vol oye (8) 15/8: 76.436 


trepilate: thaekmess (1n) 
t..; web thickness (in) 
web height (in) 

te: flange thickness (in) 
w: flange width (in) 


Vol: volume a 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

A finite element analysis program was developed and 
coupled to an optimizer for stiffened plate design. Stif- 
fened plates were designed for minimum volume with a concen- 
trated load. 

Five important conclusions are made from considering 
several design examples. 

First, a stiffened plate is much more efficient than 
the unstiffened plate for the optimum volume design. 

Second, it is required to attach the larger stiffeners 
to the diagonal elements of the plate as compared to the 
off-diagonal elements. 

Third, the T type stiffeners provide a lower volume 
design than the rectangular type stiffeners. 

Fourth, the stiffened plate with the diagonally attached 
stiffeners are much more efficient than those with the 
rectangularly attached stiffeners. 

Fifth, the stiffened plate attached diagonally with the 
levoe Stiffeners is particularly efficient for the concen- 
trated load case. 

These conclusions are based on the design of simply 
Supported plates and the conclusions could be quite differ- 


ent if other boundary conditions are considered. 


ae 





B. RECOMMENDATIONS 

The following recommendations may be useful for future 

work. 

1. The analysis program should be modified toa 
generalized program with a variety of elements. 

2. The analysis program should be extended to 
consider the dynamic problem. 

3. The analysis program should be extended to 
consider the multiple loading conditions. 

4, Routines should be added to calculate gradients 
analytically. 

5. Buckling constraints of the plate and stiffener 
elements should be added in the optimization 
process. 

6. Frequency constraints should be added in the 
optimization process. 

7. The analysis program should be modified to reduce 


the computing time during the optimization process. 
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